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ABSTRACT 


The  effects  of  the  Ionosphere  upon  radio  frequency 
electromagnetic  waves  are  discussed.  Several  types 
of  satelllte-to-satelllte  propagation  experiments  are 
presented,  and  the  problems  of  implementation  are 
delineated.  The  equipment  around  which  the  experi¬ 
ment  must  be  designed  is  described,  as  are  the  ve¬ 
hicles  available  in  Program  162.  The  relative  merits 
of  the  various  vehicles  and  the  basis  for  the  selection 
of  the  experimental  configuration  are  discussed.  The 
preliminary  design  results  are  shown. 
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Section  1 
INTRODUCTION 


Operation  of  a  satellite-to-satellite  radio  propagation  experiment  must  be  based  upon 
a  thorough  knowledge  of  the  experimental  medium,  the  experimental  environment,  and 
the  effects  thereof  upon  the  experiment.  This  report  presents  the  results  of  a  feasi¬ 
bility  study  and  preliminary  design  effort  directed  at  such  an  experiment. 

It  was  anticipated  that  the  experiment  should  include  two  transmitters,  both  in  the 
frequency  range  from  3  Me  to  20  Me,  and  would  operate  from  a  suitable  vehicle,  such 
as  the  Lockheed  Agena.  An  ejected  satellite  was  expected  to  transmit  to  the  Agena- 
type  vehicle.  The  experiment  was  anticipated  to  last  approximately  one  week. 

The  feasibility  study  was  to  consider  the  frequencies  to  be  used,  the  method  of  synchro¬ 
nizing  the  receivers  with  the  transmitters  (which  were  to  operate  alternately),  the 
pulse  length  and  repetition  frequency,  expected  characteristics  of  the  received  signals, 
relative  orbital  positions  of  the  transmitting  and  receiving  vehicles,  and  problems 
attendant  to  conducting  the  experiment  on  a  Program  162  vehicle.  The  preliminary 
design  was  to  consider  the  ejection  mechanism,  the  propulsion  rocket  (if  used),  the 
power  supply,  telemetry  system,  antennas,  beacon  transmitter,  control  unit,  and  an 
analysis  of  thermodynamic  and  other  associated  problems. 
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Section  2 

FEASIBILITY  STUDY 


2.1  IONOSPHERIC  PROPAGATION  EFFECTS* 

2.1.1  General 

An  electromagnetic  wave  propagating  through  the  atmosphere  is  affected  by  many 
phenomena.  The  variation  of  the  index  of  refraction  in  the  troposphere  (that  portion 
of  the  atmosphere  up  to  65  km  altitude)  and  in  the  ionosphere  (that  portion  of  the 
atmosphere  above  65  km  altitude)  introduces  phase  shifts  and  refraction,  causing 
apparent  deviations  in  range,  angular  position,  velocity  (Doppler  frequency),  and 
attenuation  (through  ray  path  divergence) . 

In  the  troposphere,  the  wave  is  attenuated  primarily  through  absorption  by  water  drop¬ 
lets  and  ice  formations.  In  the  ionosphere,  further  attenuation  is  produced  by  absorp¬ 
tion  and  scattering  due  to  the  appreciable  number  of  free  electrons  from  x-ray  and 
ultraviolet  radiation.  Also,  the  frequency  and  polarization  dependence  of  the  index  of 
refraction  causes  non-linear  phase  shift  effects  and  polarization  variations  (Faraday 
effect),  respectively.  Non-linear  phase  shifts  produce  distortion  (pulse  dispersion) 
in  the  information  content  of  the  electromagnetic  wave.  Because  the  scope  of  the 
satellite-to-satellite  experiment  is  concerned  mainly  with  propagation  within  the  ion¬ 
osphere,  the  effects  of  absorption,  refraction,  Faraday  rotation,  pulse  dispersion, 
and  Doppler  shifts  will  be  discussed  in  the  following  subsections. 


♦The  discussion  in  this  section  is  primarily  the  work  of  Mr.  Virgil  A.  Counter  and 
will  be  included,  in  expanded  form,  in  Space  Age  Electromagnetics  -  Volume  II: 
Propagation,  McGraw-Hill,  to  be  published  in  1963. 


2.1.2  Attenuation 


Molecular  and  Ionic  absorption  and  refractive  attenuation  are  the  primary  causes  of 
signal  attenuation  in  the  ionosphere.  The  absorption  process  involves  a  real  transfer 
of  energy  from  the  wave  via  the  electrons  to  the  gas  particles,  while  variations  in 
refractive  index  result  in  an  apparent  attenuation  of  energy. 

In  the  upper  atmosphere,  incident  x-  and  ultraviolet  radiation  are  absorbed  somewhat 
selectively  by  the  atmospheric  constituents,  resulting  in  a  more  or  less  stratified  region 
of  partially  ionized  gas,  which  is  referred  to  as  the  ionosphere.  The  lowest  region  is 
the  D- Layer  at  85  km  and  the  uppermost  is  the  r*2  Layer  at  about  290  km.  It  is  charac¬ 
teristic  of  the  ionosphere  that  its  electrical  properties  are  a  function  of  frequency.  The 
highest  frequency  which  is  reflected  by  a  given  layer  is  termed  the  critical  frequency 

u)  for  that  layer.  The  critical  frequency  is  related  to  the  local  plasma  frequency  u  , 
c  p 

which  is  defined  by 

ojp  =  56.3  VN  (inks  units)  (2.1) 

where  N  =  local  electron  density. 


The  relation  between  u>  and  u>  is  given  by 

c  P 


=  w  =  56.3./T? 

pm  ax  V  max 


(2.1a) 


where  the  subscript  max  refers  to  the  maximum  for  the  layer  in  question.  The  reflec¬ 
tion  process  is  very  similar  to  that  found  in  metals.  The  ionospheric  electrons  are 
caused  to  move  by  the  incident  field,  but  their  motion  lags  the  electric  field  by  90  deg 
(in  the  lossless  case),  while  the  displacement  current  leads  by  90  deg.  When  the  two 
currents  are  equal  in  magnitude,  reflection  is  complete.  If,  however,  there  are 
numerous  neutral  and  ionized  gas  particles  (as  indeed  there  are  in  the  ionosphere) , 
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then  a  second  process  is  involved.  The  electrons  will  now  suffer  a  number  of  collisions 
each  second  and,  with  each  collision,  the  electrons  will  lose  part  of  its  kinetic  energy. 
Thus,  while  reflection  will  still  occur,  the  process  is  no  longer  lossless  so  far  as  the 
electromagnetic  field  is  concerned. 

Refractive  attenuation  results  from  the  divergence  of  the  ray  paths,  which  is  caused  by 
the  variation  of  the  index  of  refraction.  The  gross  structure  of  the  atmosphere  con¬ 
tributes  to  that  portion  of  the  refractive  attenuation  which  usually  is  positive  and  actually 
attenuates  the  signal.  However,  for  situations  in  which  the  ionosphere  is  traversed  by 
the  signal,  the  refractive  attenuation  can  be  negative  and  the  signal  can  be  enhanced  by 
the  effect  of  refraction.  Local  inhomogeneities  can  cause  this  refractive  attenuation  to 
fluctuate  positively  and  negatively  at  amplitudes  of  many  db  at  VHF  and  lower  frequen¬ 
cies.  Therefore,  at  certain  periods  during  the  fluctuations,  signal  enhancement  caused 
by  refraction  may  more  than  offset  the  attenuation  due  to  absorption  and  scattering. 

As  implied  in  the  above  discussion  on  absorption,  the  critical  frequency  and,  therefore, 
the  refractive  ionospheric  attenuation  are  approximately  proportional  to  ihe  inverse 
square  of  the  frequency  and  generally  are  negligible  above  1000  me  for  the  most  severe 
ionospheric  conditions.  At  a  frequency  of  100  Me,  attenuation  due  to  absorption  is  less 
than  1  db,  but  refractive  fluctuations  may  be  greater  than  10  db.  As  the  radio  frequency 
approaches  the  ionospheric  critical  frequency,  both  the  absorptive  and  refractive  space- 
to-ground  attenuations  may  become  tens  of  db. 

It  has  been  shown*  that  the  attenuation  constant,  a,  is  given  by 


(nepers/m) 

where  v  is  the  electron  collision  frequency  and  C  is  the  velocity  of  light. 

*E.  C.  Jordan,  Electromagnetic  Waves  and  Radiating  Systems,Prentice  -  Hall,  Inc. 
p.  673.  -  - - 
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Figure  2-1  presents  this  relationship  in  a  normalized  form  with  v/u  versus  u>/w 

P  P 

with  a/w  as  a  parameter.  As  an  example,  consider  the  case  where  f  -  16  Me 

(w  =  108),  f  =4Mc(w^2.5  x  10^),  v  =  1.5  x  10**.  For  these  conditions  it  is  found 
r  -14  -6  —3 

that  a/io  is  roughly  6  x  10  and  a  =  1.5  x  10  db/m,  or  2.8  x  10  db/nm. 

P  *7 

The  effect  of  reducing  frequency  is  seen  if  f  is  changed  to  8  Me  (w  =  5  x  10  );  this 


yields  a 


10  ®  db/m  or  1.85  x  10  ^  db/nm. 


For  the  low  collision  frequencies  typical 


of  the  ionosphere  (see  Fig.  2-2),  a  rises  quite  abruptly  as  w  approaches  from 
above. 


2.1.3  Faraday  Effect 

The  earth's  magnetic  field  causes  the  electrons  of  the  ionosphere  to  deviate  from  their 
normal  oscillatory  paths  when  excited  by  the  electric  field  of  an  electromagnetic  wave. 
Therefore,  the  polarization  of  the  forward  scattered  field  from  each  electron  Is  shifted 
slightly  from  that  of  the  incident  field.  This  accumulative  shift  of  polarization  as  the 
wave  propagates  through  the  Ionosphere  is  called  the  Faraday  effect. 

One  aspect  of  this  effect  is  the  longitudinal  Faraday  effect,  or  Faraday  rotation,  which 
is  a  rotation  of  the  polarization  ellipse  of  an  arbitrarily  polarized  plane  wave.  The 
amount  of  rotation  depends  on  the  component  of  the  magnetic  field  along  the  wave  normal 
and  varies  approximately  as  the  inverse  square  of  the  wave  frequency.  At  frequencies 
below  200  Me,  the  transverse  Faraday  effect  also  is  significant,  but  only  when  the  wave 
normal  is  close  to  being  perpendicular  to  the  direction  of  the  magnetic  field.  This  effect 
can  change  the  eccentricity  of  the  polarization  ellipse,  causing  a  linearly  polarized  wave 
to  become  circularly  polarized  (or  vice  versa),  or  causing  a  circularly  polarized  wave 
to  alter  its  sense.  The  transverse  effect  is  proportional  to  that  component  of  the  mag¬ 
netic  field  perpendicular  to  the  wave  normal  and  varies  approximately  as  the  inverse  cube 
of  the  frequency. 

A  common  method  of  analysis  of  the  Faraday  effect  is  to  consider  the  original  wave  as 
splitting  into  two  magneto-ionic  components  designated  as  the  ordinary  and  extra¬ 
ordinary  waves.  The  two  components  propagate  independently  with  different  phase 
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velocities.  The  polarizations  of  the  two  components  are  orthogonal  when  absorption  is 
not  a  problem.  If  the  original  signal  is  modulated  and  the  delay  difference  between  the 
two  components  is  large  compared  to  the  periods  of  modulation,  the  Faraday  effect 
would  cause  time-delay  distortion  of  the  received  signal  (unless  only  a  single  component 
is  monitored).  In  pulsed  transmissions,  this  distortion  would  increase  the  length  or 
cause  multiplicity  of  the  pulses. 

The  magnitude  of  the  Faraday  effect  may  be  calculated  in  terms  of  thp  relative  phase 
shift  between  the  two  modes,  given  by 

A^F  =  f  /  (An>  ds  (2-3) 


where  An  =  n  -  n„  =  the  difference  in  the  indices  of  refraction  for  the  two  modes  of 
o  E 

propagation.  This  quantity  may  be  written  as 


An  =  n  - n 
o  Jb 


when  n  and  n_  differ  considerably,  or,  as 
o  r. 


A  n 


„  2/2 
2u>  /w 

=  E L 

n„  +  n_ 
o  £ 


- 

2 

- 

■ 

‘4/4  .  4  ,  2 

sin  ^  WH  2 

(wr/w)  sin2  4> 

2 

WH  2 

,  r,  +  cos  ib 

i  _  U) 

1  "  2  2 

-  2  008  * 

/  2  2\ 

/w  -wl 

U)  -  O) 

P 

u> 

2 

L  CO 

(2 : 

when  nQ  w  n^  and  where  is  the  gyro  frequency  and  $  is  the  angle  of  the  wave 
normal  to  the  magnetic  field.  For  the  quasilongitudinal  case,  when  nQ  and  nE  are 
very  close  to  1,  and  ^  is  approximately  zero 


A< p 


-1  /  “p  “H  \ 

c  V2  -  4J 


ds  rad/m 


(2.5) 
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and  when  w  and  u>„  are  constant 
P  H 


A 


OJ  +  U) 


H 


2  2 
w  -  wH 


R 


(2.  5a) 


where  R  is  the  path  length. 

In  the  quasi-transverse  case,  when  n  and  n,,  are  very  close  to  1 ,  and  4  is 

o  hj  T 

approximately  r/2 


A  <p 


(2.6) 


and  when  w 

P 


and  u 


H 


are  constant 


A  (p 


_1 _ 

2  w  C 


2  2 

U) 

P  « 


R 


u) 


Gu) 


U> 


H 


(2.  6a) 


Taking  the  case  previously  discussed  (f  =  16  Me,  f  =4  Me)  and  using  f„  =  1.3  Me, 

P  h 

(40  deg  north  latitude),  the  magnitude  of  the  longitudinal  effect  is  1.7  x  10"3  rad/m 

(0. 1  deg/m)  and  the  magnitude  of  the  transverse  effect  is  1. 5  x  10-4  rad/m  (0. 009  deg/m) 

As  it)  approaches  the  magnitude  of  these  effects  increases  rapidly. 

The  magnitude  of  the  rotation  of  the  polarization  vector  in  the  quasi- longitudinal  case 
is  one-half  the  phase  shift  between  the  two  components.  In  the  quasi-transverse  case, 
however,  the  effect  is  rather  complex  and  depends  upon  the  polarization  of  the  trans¬ 
mitted  wave  relative  to  the  magnetic  field  of  the  earth. 
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2.1.4  Pulse  Dispersion 


If  the  phase  delay  through  a  medium  increases  linearly  with  frequency,  the  signal 
modulation  is  only  shifted  in  time,  with  no  distortion.  Therefore,  any  medium  which 
introduces  non-linear  phase  delays  will  distort  any  modulation  on  a  passing  signal. 
This  distortion  may  be  estimated  with  the  aid  of  an  analysis  by  Elliott*  and  some 
reasonable  generalizations.  Elliott  derived  an  equation  for  the  shape  of  a  single, 
originally  rectangular  pulse  after  it  has  traveled  through  a  lossless  section  of  wave¬ 
guide.  He  assumed  that  the  phase  constant  is  proportional  to  the  phase  delay  and  can 
be  represented  by  the  first  three  terms  of  the  Taylor  series.  With  these  conditions, 
Elliott's  expression  can  apply  to  the  distortion  of  a  rectangular  pulse  in  any  dispersive 


homogeneous  medium  for  which 


<  0  and  wT  »  1 . 


His  results  were  expressed  as 


y 


(2.7) 


where 


4>  =  actual  phase  delay  over  the  path  of  wave  propagation 

f  =  carrier  frequency 
u>  =  angular  carrier  frequency 
T  =  time  duration  of  the  pulse 
c  =  velocity  of  light  in  free  space 

Q  =  80.6  JNdR=  a  function  proportional  to  the  integrated  electron  density 
along  the  path  of  wave  propagation 

N  =  electrons  per  cubic  meter 

dR  =  meters 


*R.  S.  Elliott,  "Pulse  Waveform  Degradation  due  to  Dispersion  in  Waveguide,"  IRE 
Transactions  on  Microwave  Theory  and  Techniques,  Vol.  MTT-5,  No.  4.  Oct  1957, 
pp.  254-257. 
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Thus,  for  a  given  pulse  width,  frequency  and  Q ,  the  amount  of  pulse  distortion  may  he 
approximated  by  using  the  above  equation.  Figure  2-3  shows  the  pulse  Bhape  for  several 
values  of  y . 

Elliott's  results  also  may  be  applied  qualitatively  to  a  train  of  pulses.  It  is  estimated 
from  Fig.  2-3  that,  for  y  =  1 ,  the  effective  signal-to-noise  ratio  for  a  pulse  train 
would  be  reduced  by  approximately  3  db.  In  fact,  the  pulse  degradation  with  this  value 
of  y  indicates  the  difficulty  in  the  resolution  of  pulses  with  a  time  between  pulses 
equal  to  the  pulse  widths.  Therefore,  it  appears  that  the  maximum  distortion  permis¬ 
sible  for  a  high  data  rate  pulse  modulation  is  when  y  =  1 .  In  the  case  of  an  inhomo¬ 
geneous  and  turbulent  ionosphere,  this  value  may  be  much  less. 

Considering  the  specific  case  where  f  =  16  Me,  T  =  250  hb  ,  prf  =  10 ,  and  the 
electron  density  is  constant  at  2  x  10*  */M^  (fp  =  4  Me),  it  is  found  that  y  =  1.27  x  10”® 
Therefore,  the  criterion  cited  above,  y  =  1 ,  will  be  obtained  when  the  separation  reachei 
79  km  (42. 5  nm )  and  it  must  be  expected  that  at  the  peak  of  the  F2  region,  pulse  shape 
degradation  will  cause  a  serious  drop  in  the  instantaneous  peak  power  received. 

2.1.5  Doppler  Effect 

The  relative  radial  velocity  component  of  a  transmitting  or  receiving  source  can  be 
determined  from  the  frequency  difference  (or  Doppler  frequency)  shift  between  the  two 
signals.  In  free  space,  this  Doppler  frequency  can  be  expressed  as 


where 

fr  =  received  frequency 

fD{8  *  Doppler  frequency  in  free  space 

R  =  relative  radial  velocity  between  the  transmitter  and  receiver 
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During  actual  measurements,  this  equation  would  represent  the  average  one-way 
frequency  difference  over  a  finite  time  interval. 

For  the  dispersive  ionosphere,  the  apparent  relative  radial  velocity  will  be  a  function 
of  frequency  and  the  Doppler  equation  must  be  written 


if  the  collision  frequency  is  small  relative  to  f .  Use  of  two  or  more  frequencies, 
coherent  with  one  another,  will  permit  measurement  of  the  Doppler  component  result¬ 
ing  from  the  dispersive  nature  of  the  medium,  and  thus  will  yield  the  Integrated  elec¬ 
tron  density  between  the  transmitters  and  the  receivers.  In  the  case  where  N  is 
constant  along  the  path,  the  Doppler  frequency  becomes 


(2.9a) 


where  it  must  be  remembered  that  n  contains  a  factor  of  4  and  then  the  frequency 
dependence  and  range  change  rate  dependence  are  obvious. 


An  idea  of  the  magnitude  of  the  dispersive  Doppler  effect  for  slowly  separating  satel¬ 
lites  can  be  obtained  from  the  following  calculation.  Assume  that  f ^  =  16  Me , 

f  =  8Mc,  f  =  4  Me ,  and  R  =  8  m/sec. 

I  p 

f  =-4.3  cps  at  16  Me ,  -  2. 16  cps  at  8  Me 
Dfs 

f  =-4.15  cps  at  16  Me ,  -  1. 86  cps  at  8  Me 
F 

m 
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The  Doppler  components  due  to  the  medium  are  +  0. 15  cps  at  16  Me  and  +  0.  29  cps 
at  8  Me  -  certainly  not  very  large  changes. 

Thus,  while  the  dispersive  Doppler  technique  has  been  highly  successful  for  ground 
to  rocket  measurements  of  integrated  electron  density,  it  is  not  a  very  good  tool  for 
measuring  integrated  electron  density  between  two  slowly  separating  satellites.  This 
is  not  to  imply  that  satellite-to-ground  measurements  are  not  feasible;  these  measure¬ 
ments  can  be  performed  much  as  in  the  case  of  the  vertical  rocket.  The  data  reduction 
problem  is  more  complex,  however,  because  the  data  received  relates  to  a  number  of 
ray  paths,  not  just  one.  Furthermore,  when  the  transmitting  frequencies  are  close  to 
local  plasma  frequency,  the  ray  paths  will  be  sufficiently  different  to  degrade  the  data 
quite  seriously. 

2.2  EXPERIMENTAL  APPROACHES 

2.2.1  Sounding  Approach 

Previous  radio  frequency  studies  of  the  ionosphere  have,  of  necessity,  been  confined 
either  to  a  ground-to-ground,  or  a  space-to-ground  propagation  mode  of  operation. 
Despite  the  limitations  imposed  by  these  restrictions,  a  number  of  effective  experi¬ 
mental  techniques  have  been  developed  and  a  wealth  of  data  has  been  accumulated. 

A  large  fraction  of  the  experimental  ionospheric  data  has  been  acquired  from  vertical- 
or  oblique-incidence,  swept-frequency  pulse  soundings.  Such  techniques  are  equally 
applicable  for  "top  sounding"  from  a  satellite  (or  satellites)  and  can  prove  of  great 
value  in  describing  the  upper  reaches  of  the  ionosphere. 

2.2.2  Phenomenological  Approach 

A  greater  degree  of  freedom  is  afforded  by  transmission  from  satellite  to  satellite; 
this  makes  possible  an  entirely  new  type  of  operation,  wherein  both  the  transmitter 
and  the  receiver  are  located  within  the  ionosphere.  Under  these  conditions,  it  is 
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possible  to  study  the  ionospheric  propagation  effects  mentioned  in  Section  2. 1  above. 

The  experiment  can  yield  more  useful  data  if  the  satellites  are  placed  in  an  orbit  such 
that  the  vehicles  rise  and  fall  in  unison  through  the  ionosphere.  Because  of  the  general 
incompatabiiity  of  the  top-sounder  and  the  propagation  effect  approaches,  they  will,  in 
this  study,  be  considered  separately;  and  though  regions  of  overlap  may  exist,  they 
generally  will  be  ignored  in  this  report. 

2.3  EXPERIMENTAL  TECHNIQUES 

2.3.1  Experimental  Geometries 

A  wide  variety  of  techniques  may  be  used  to  implement  the  two  general  approaches 
discussed  above;  for  the  sake  of  clarity,  description  of  these  techniques  will  also  be 
kept  quite  general  in  the  following. 

A  top-sounder  experiment  may  be  designed  as  a  one-satellite  operation,  but,  in  the 
case  where  the  experiment  is  cf  secondary  importance  to  the  satellite  flight,  it  may 
prove  necessary  to  separate  the  transmitter  from  the  receiving  satellite  to  avoid 
interference  with  primary  payload  operation.  Such  a  two-vehicle  experiment  becomes, 
of  necessity,  an  oblique- incidence  sounder.  A  sounding  experiment  should  be  conducted 
from  vehicles  well  above  the  peak  of  the  ionosphere;  thus  an  altitude  of  at  least  500  to 
1,000  nm  is  indicated. 

A  two-satellite  ionospheric  phenomena  experiment  may  take  several  forms,  depending 
upon  the  phenomena  to  be  studied.  Similar  to  the  top  sounder,  but  at  ionospheric  alti¬ 
tudes,  an  experiment  could  consist  of  two  satellites  relatively  close  together.  Such  a 
system  could  obtain  data  on  Faraday  rotation,  multipath  effects,  pulse  dispersion, 
attenuation,  and,  to  a  minor  degree,  ducting.  In  contrast,  use  of  antipodal  satellites 
would  emphasize  study  of  long  range  ducting  effects  and  could  obtain  little  or  no  data  on 
other  propagation  phenomena. 
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Between  these  extremes  any  number  of  orbital  geometries  can  be  devised,  each  with 
its  relative  merits.  Worth  consideration  is  the  set  of  geometries  in  which  the  vehicles 
vary  in  separation  from  minimum  to  maximum  and  vice  versa  a  number  of  times  during 
the  life  of  the  experiment.  The  number  of  times,  of  course,  may  be  adjusted  as  desired. 

2.3.2  Vehicular  Considerations 

Not  only  must  a  variety  of  geometries  be  considered  in  the  experimental  design,  but 
the  various  ways  in  which  these  geometries  can  be  achieved  must  be  studied.  This 
discussion,  also,  will  be  of  a  general  nature  and  will  not  deal  with  specific  vehicles. 

It  is  technically  quite  difficult  to  conduct  an  oblique- incidence  top-sounding  experiment 
between  two  satellites  which  are  not  placed  in  orbit  by  the  same  vehicle.  However,  if 
the  satellite  orbits  do  derive  from  the  same  launch,  then  vehicle  separation  can  be 
achieved  only  by  means  of  a  continuous  change  in  separation  as  the  experiment  progresses. 
Use  of  two  launch  vehicles  to  achieve  a  more  or  less  fixed  separation  between  satellites 
may  be  considered;  but  such  a  feat  is  roughly  as  difficult  as  achieving  a  rendezvous  and 
requires  precise  synchronization  on  the  ground  and  some  degree  of  station-keeping 
capability  in  one  of  the  satellites.  It  should  be  noted  that  a  pure  antipodal  experiment 
requires  almost  as  precise  ground  synchronization  as  the  stabilized  oblique  sounder. 

As  mentioned  previously,  the  experiment  may  be  designed  with  a  variable  geometry;  and 
this  mode  of  operation  can  be  obtained  in  several  ways.  The  two  satellite  orbits  may 
derive  from  a  single  launch  and  their  initial  separation  rates  may  be  increased  from 
several  feet  per  second  to  several  hundred  feet  per  second.  Thus  the  time  required  to 
achieve  vehicle  opposition  would  range  from  several  months  to  several  days,  depending 
upon  altitude,  with  the  lowest  velocities  providing  a  quasi-stationary  state  of  operation. 

The  separate-launch,  variable-geometry  approach,  where  little  if  any  attempt  is  made 
to  control  orbital  parameters,  will  in  general  provide  for  an  infinity  of  geometries. 

This  approach  also  admits  of  multiple  vehicle  operation,  including  the  special  case 
wherein  one  vehicle  has  a  sufficiently  long  operating  life  to  operate  with  a  number  of 
short-lived  vehicles  launched  sequentially.  Such  a  possibility  will  be  discussed  at 
greater  length  in  Section  3.3.2. 
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2.4  CONCLUSIONS 


It  has  been  shown  that  the  propagation  of  high-frequency  signals  from  one  satellite  to 
another  will  be  affected  in  several  ways.  These  effects  are  of  sufficient  maghitude  as 
to  be  observable  even  when  using  relatively  simple  equipment  on  the  satellites.  Not 
only  may  ionospheric  propagation  phenomena  be  observed  when  transmitting  from  satel¬ 
lite  to  satellite,  but  it  is  also  possible  to  operate  a  two-satellite  ionospheric  top  sounder. 

It  should  be  understood  that  no  experimental  geometry  permits  studying  all  aspects 
of  the  ionosphere  at  one  time.  With  some  variable  geometries,  however,  most  iono¬ 
spheric  phenomena  can  be  studied  at  one  time  or  another  during  the  lifetimes  of  the 
vehicles.  Although  almost  any  experimental  geometry  desired  may  be  obtained,  some 
are  much  easier  to  achieve  than  others  and  this  may  prove  to  be  a  controlling  factor  in 
the  design  of  the  experiment.  Thus,  it  may  be  concluded  that  any  one  of  a  number  of 
ionospheric  propagation  experiments  are  feasible;  but  prior  to  the  selection  of  the 
experiment  to  be  undertaken,  it  is  necessary  to  weigh  with  great  care  the  relative  merits 
of  both  the  possible  experiments  and  the  methods  available  for  implementation  of  the 
experiments.  The  actual  selection  was  part  of  the  preliminary  design  phase  and  is 
discussed  in  the  following  section. 
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Section  3 

PRELIMINARY  DESIGN 


3.1  GENERAL 

The  preliminary  design  of  the  satellite-to-satellite  radio  propagation  experiment  will 
be  controlled  by  certain  initial  constraints.  It  was  predetermined  that  two  frequencies 
in  the  3  to  20  Me  region  should  be  used.  The  experiment  will  use  a  suitable  vehicle  to 
launch  the  second  vehicle  into  an  orbit  approximately  in  the  same  plane  as  the  mother 
vehicle,  with  the  transmitter  in  the  ejected  package.  It  is  anticipated  that  the  mother 
vehicle  will  be  an  Agena  and  that  the  experimental  lifetime  will  be  about  one  week.  As 
part  of  the  second  item  of  work  (design  and  construction  of  an  operating  model)  it  is 
specified  that  the  Government  should  furnish  the  transmitters,  receivers,  power  con¬ 
verters,  and  antenna  tuning  units;  these  items,  designed  and  built  by  General  Electric 
Company,  turther  define  the  experimental  design.  Several  possible  launching  schemes 
will  be  investigated  because  the  original  concept  appears  to  be  too  restrictive.  Other 
than  this  one  deviation,  the  preliminary  design  will  proceed  roughly  as  originally 
planned. 

3.  2  GOVERNMENT-FURNISHED  EQUIPMENT 

The  Government-furnished  General  Electric  transmitting  and  receiving  units  are  solid- 
state,  high-efficiency  units  designed  to  meet  the  temperature,  vibration,  and  shock  en¬ 
vironments  encountered  in  missile  and  satellite  flights.  The  system  consists  of  a  trans¬ 
mitter  programmer,  two  transmitters,  two  receivers,  and  an  IF  amplifier.  The  trans¬ 
mitters  are  frequency-controlled  by  a  crystal  and  have  a  peak  powa:  output  of  250  w  up  to 
10  Me  (150  w,  from  10  to  20  Me).  The  nominal  pulse  width  is  250  psec,  the  prf  is  10, 
and  the  average  peak-power  requirement  of  the  transmitter  is  1. 9  w.  The  transmitter 
programmer  unit  consists  of  stabilized  digital  count-down  circuits,  which  generate  the 
correct  prf  and  pulse  width,  and  timing  circuits,  which  select  the  correct  transmitter. 
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Its  power  requirement  is  0.72  w  The  receivers  are  single-conversion  superheterodyne 
units  with  the  gain  and  noise  figure  commensurate  with  the  received  noise  level.  The 
power  requirement  of  the  receiver  is  0.  17  w,  and  the  output  frequency  is  3  Me.  The  IF 
amplifier  has  a  10-ke  bandwidth  centered  at  3  Me,  a  minimum  gain  of  72  db,  an  age  range 
of  50  db,  and  age  and  video  outputs.  The  power  requirement  of  the  IF  strip  is  0.44  w. 

All  units  require  a  28-v  power  source  except  the  programmer,  which  requires  -12  v 
and  +16  v,  and  the  transmitter,  which  requires  +16  v,  +28  v,  and  +50  v.  The  "non¬ 
standard"  voltages  will  be  obtained  from  a  28-v  power  converter;  however,  this  unit  will 
not  be  GFE  as  required  by  the  original  contract. 

It  is  anticipated  that  both  receiver  and  transmitter  specifications  will  be  changed  some¬ 
what  prior  to  delivery  of  these  items,  but  the  above  represents  presently  available  data 
furnished  and  confirmed  by  the  manufacturer. 

3.  3  EXPERIMENTAL  CONFIGURATIONS  OBTAINABLE 

Because  of  the  dependence  of  experimental  geometry  upon  specific  launch  vehicle  details, 
the  preliminary  design  effort  was  predicated  upon  the  use  of  the  basic  Agena  vehicle, 
as  used  in  the  LMSC  Program  162.  It  is,  therefore,  essential  that  operations  of  this 
program  be  studied.  It  is  the  practice  of  the  Air  Force  to  utilize  excess  space  and 
weight  allowance  on  Program  162  vehicles  for  research  payloads.  These  secondary 
payloads  are  placed  in  the  Agena  vehicle  and  carried  into  orbit  on  a  non-interference 
basis.  Research  payload  data  are  generally  fed  via  two  sixty-point-per-second  com¬ 
mutators  to  a  two-channel  tape  recorder.  Upon  ground  command,  the  tape-recorded 
data  is  played  back  at  high  speed  through  the  telemetry  link  It  is  possible,  following 
completion  of  the  primary  mission,  to  continue  operation  of  the  research  payload  with 
the  primary  vehicle  functions  shut  down.  By  this  means  the  experiment  may  be  pro¬ 
longed  for  a  period  of  a  week  or  so.  Such  a  vehicle  is  frequently  referred  to  as  a 
"zombie."  A  further  aspect  of  Program  162  is  the  planned  installation  of  a  large, 
self-propelled  satellite  on  the  Agena.  This  satellite,  referred  to  as  the  P-11,  is  to 
be  ejected  once  the  Agena  is  in  orbit,  and  it  will  then  propel  itself  into  one  of  two 
higher  orbits.  The  typical  Program  162  Agena  has  a  perigee  of  115  nm  and  an  apogee 
between  225  nm  and  380  nm  and  the  P-11  will  be  in  either  a  200  nm  circular  orbit 
or  in  an  orbit  with  a  perigee  of  115  nm  and  an  apogee  of  2, 500  nm. 
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On  the  basis  of  these  operations  it  is  possible  to  define  three  basic  modes  of  operation 
for  a  satellite-to-satellite  experiment:  (1)  two  Agenas,  (2)  an  Agena  and  a  P-11,  or 
(3)  an  Agena  and  an  ejected  package.  (The  possibility  of  ejecting  a  package  of  the  size 
envisioned  here  from  a  P-11  must  be  rejected  because  the  general  construction  of  the 
P-11  precludes  such  an  operation.) 

3.3.1  Concept  1 

This  concept,  using  two  Agena  "zombies,"  Involves  a  rather  complex  and  unpredictable 
geometry.  Vehicles  in  Program  162  are  placed  into  slightly  elliptical  orbits  and  the 
apogee  is  intended  to  be  about  225  nm  but  occasionally  it  has  been  as  high  as  380  nm. 
This  variation  in  apogee  is  random  in  nature,  as  is  a  variation  in  orbital  period.  The 
time  of  vehicle  launch  also  is  subject  to  large  variations  due  to  unforeseen  problems 
encountered  during  launch  operations.  For  these  two  reasons,  it  is  not  possible  to 
predict  in  advance  the  inter-satellite  geometry  to  be  expected  during  the  course  of 
the  experiment.  It  must  be  expected,  for  this  concept,  that  during  the  experimental 
period  the  separation  would  be  as  great  as  one-half  the  earth's  circumference  and 
might  be  as  small  as  a  few  hundred  miles  and  probably  would  vary  within  these  limits 
several  times  during  the  life  of  the  vehicles 

Given  this  experimental  geometry,  it  is  necessary  that  the  experiment  function  at 
separations  of  up  to  11, 000  nm.  Postulating  the  experimental  study  objective  to  be 
ionospheric  ducting,  multi-path  communications  and  interference,  antipodal  effects, 
etc. ,  it  appears  advisable  to  set  one  frequency  several  times  greater  than  the  maxi¬ 
mum  anticipated  critical  frequency  of  the  F2  region  of  the  ionosphere  to  preclude 
Incurring  excessive  signal  loss.  To  study  the  same  effects  at  times  of  significantly 
reduced  ionospheric  activity,  the  second  frequency  should  be  at  least  half  the  first 
frequency  and  yet  should  still  be  above  the  critical  frequency  of  the  F2  region  most 
of  the  time.  To  this  end,  it  is  suggested  that  the  frequencies  be  roughly  7  Me  and 
14  Me.  The  transmitting  package  would  consist  of  a  pair  of  antennas,  an  antenna 
tuning  unit,  two  transmitters,  a  programmer,  a  power  converter,  interconnecting 
wires,  a  container,  and  a  control  system  which  will  activate  the  antenna  erection 
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mechanism  and  the  programmer.  The  transmitters  would  be  keyed  by  the  transmitter 
programmer  and  would  transmit  250-psec  pulses  with  a  peak  power  of  about  150  w 
and  a  prf  of  10.  Kaeh  transmitter  would  operate  for  roughly  1/2  sec  followed  by  1/2 
see  of  no  transmission,  after  which  the  other  transmitter  would  be  keyed.  The  keying 
could  be  coded  for  frequency  identification.  The  power  converter  would  be  used  to 
convert  vehicle  28  vdc  to  the  dc  voltages  required  for  the  transmitters  and  the  pro¬ 
grammer.  Total  power  consumption  for  the  transmitting  package  would  be  2.  6  w. 

The  receiving  package  would  consist  of  a  pair  of  antennas,  two  receivers,  an  IF  ampli¬ 
fier,  interconnecting  wires,  a  container,  and  a  control  mechanism  which  would  turn 
on  power  to  the  receivers  and  the  IF  amplifier.  Both  receivers  feed  the  IF  amplifier 
simultaneously;  therefore,  it  would  be  necessary  to  turn  off  the  receiver  not  receiving 
the  desired  signal.  The  General  Electric  receivers  would  require  a  programmer  to 
accomplish  this.  All  receiving  units  require  28  vdc,  and  the  power  consumption  of  the 
package  will  be  0.75  w.  The  receiver-IF  amplifier  outputs  are  video,  AGO,  and  peak 
envelope,  and  it  would  be  planned  to  telemeter  video  to  the  Air  Force  ground  tracking 
sites  on  a  real-time  basis  during  the  command  telemetry  period  and  to  record  AGC  and 
peak  envelope  for  command  dump. 

Total  weight  for  the  transmitter  is  about  20  lb  and  for  the  receiver  about  15  lb.  Inasmuch 
as  both  vehicles  would  be  tracked  by  Air  Force  telemetry  stations,  all  information  neces¬ 
sary  to  determine  the  inter-satellite  geometry  would  be  available.  These  data,  the 
telemetered  data,  and  ionospheric  activity  data  derived  from  vertical  soundings  and 
other  sources  would  be  used  to  study  the  various  aspects  of  ionospheric  propagation. 

The  experiment  would  also  provide  signals  which  might  be  used  by  other  scientific  or¬ 
ganizations  to  study  the  many  space-to-ground  propagation  phenomena  which  have  been 
the  subject  of  considerable  theoretical  study  but  little  experimental  work. 

3.3.2  Concept  2 

The  P-11  experiment,  unlike  the  "zombie"  concept  discussed  above,  would  be  reasonably 
well  defined  in  advance  insofar  as  inter-satellite  geometry  is  concerned.  The  Agena 
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orbit  would  vary  as  previously  noted,  and  the  P  -11  would  enter  a  circular  orbit  at  an 
altitude  of  about  200  nm;  the  two  orbits  would  be  co-planer.  Because  of  the  different 
periods  to  the  orbits,  the  two  vehicles  would  separate  quite  rapidly  and  soon  be  separated 
by  11,000  nm.  Upon  reaching  maximum  separation,  the  range  between  the  satellites 
would  close,  and  the  cycle  would  then  begin  anew.  It  is  readily  understood  that  this 
geometry  is  considerably  less  complex  than  that  involved  in  the  "zombie"  experiment. 

The  transmitting  P-11  would  last  at  least  six  months  and  thus  could  be  used  by  a  number 
of  receiving  satellites,  allowing  considerable  flexibility  of  operation. 

The  extreme  range  involved  requires  that  the  experiment  resemble  the  "zombie"  ex¬ 
periment  in  many  respects.  The  fact  that  the  P-11  would  remain  above  the  peak  of 
the  F2  layer  makes  it  desirable  to  have  the  lower  frequency  somewhat  lower  than  that 
proposed  for  the  "zombie"  so  that  some  details  of  the  upper  "side"  of  the  ionosphere 
might  be  studied  when  the  vehicles  are  close  together.  To  perform  such  measure¬ 
ments  the  lower  frequency  should  be  close  to  the  critical  frequency  of  the  ionosphere; 

4  to  5  Me  would  probably  be  the  frequency  range  involved. 

The  contents  of  the  transmitting  and  receiving  packages  would  be  essentially  the  same 
as  those  previously  described,  except  for  changes  due  to  the  different  low  frequency 
to  be  used  in  the  P-11  configuration  and  such  changes  in  sequencing  as  are  deemed 
necessary.  The  transmitting  antennas  might  have  to  be  separate  from  the  transmitter 
package,  and  the  power  for  the  transmitters  would  be  derived  from  the  P-ll's  solar 
cells;  but  operation  of  the  equipment,  recovery  of  the  data,  and  the  weights  would  be 
as  proposed  for  the  "zombie"  experiment. 

3.3.3  Concept  3 

This  concept  involves  a  sub-satellite  designed  to  separate  from  a  Program  162  vehicle, 
carrying  with  it  the  transmitting  package.  The  sub-satellite  would  contain  its  own 
power  supply  and  a  beacon  for  tracking  purposes,  and  for  technical  reasons,  would  sepa¬ 
rate  from  the  Agena  at  a  very  low  velocity.  The  mother  vehicle  would  have  the  same 
orbit  discussed  previously  and  because  of  the  low  separation  velocity,  the  sub-satellite 
would  have  essentially  the  same  orbit  as  the  mother  vehicle. 
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Because  of  the  orbit  eccentricity,  the  two  satellites  would  operate  above,  in,  and  below 
the  maximum  of  the  F9  layer  at  rather  short  ranges.  It  is,  therefore,  suggested  that 
the  experiment  should  study  ducting,  ionospheric  attenuation,  pulse-shape  distortion 
(frequency  dispersion),  ionospheric  reflection,  etc.  For  those  purposes  it  appears 
test  to  operate  at  frequencies  of  7  and  14  Me. 

The  transmitting  sub-satellite  would  include  the  components  discussed,  plus  batteries 
and  ejection-spin  mechanism.  The  receiving  package  would  be  the  same  as  in  the  other 
concepts.  The  transmitting  satellite  will  weight  at  least  50  lb  and  the  receiving  package 
will  weigh  about  15  lb. 

Only  the  Agena  would  be  tracked  by  the  Air  Force  telemetry  stations;  the  services  of 
SPADETS  and  similar  organizations  would  be  enlisted  to  track  the  ejected  satellite. 
Acquisition  of  the  experimental  data  would  be  the  same  as  the  other  concepts,  and 
auxiliary  scientific  studies  could  be  performed  as  previously  described. 

3.4  EXPERIMENT  SELECTION 

It  is  necessary,  as  part  of  the  preliminary  design,  that  the  type  of  experiment  to  be 
performed  be  determined.  The  preceding  section  has  discussed  the  experiments  pos¬ 
sible  within  the  Program  162  framework.  The  relative  advantages  and  disadvantages 
of  these  approaches  arc  presented  below  for  comparison. 

3.4.1  Zombie  to  Zombie 

Advantages: 

(1)  The  weight  penalty  per  Agena  is  a  minimum  for  this  concept. 

(2)  The  possibility  of  obtaining  vehicle  space  is  quite  high. 

(3)  The  costs  of  the  experimental  packages  are  about  the  same  as  for  P-11 
and  less  than  the  sub- satellite. 

(4)  A  wide  range  of  inter-satellite  geometries  can  be  anticipated. 

(5)  The  experimental  life  would  be  roughly  three  weeks. 
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Disadvantages: 

(1)  The  experimental  geometry  is  very  complex  and  data  reduction  will  thus  be 
complicated. 

(2)  The  probability  of  experiment  success  is  possibly  slightly  lower  than  for 
the  other  concepts. 

(3)  The  radiation  from  the  transmitters  might  present  a  very  serious  RFI 
problem  for  the  other  experiments  on  board. 

(4)  There  is  no  way  at  present  to  assure  that  the  two  Agenas  would  be  launched 
sequentially;  therefore,  it  is  possible  that  both  launches  would  be  successful 
but  so  far  apart  in  time  as  to  preclude  cooperation  between  vehicles. 

3.4.2  P-11  to  Agena 

Advantages: 

(1)  This  concept  is  only  slightly  heavier  in  total  weight  than  the  Zombie 
experiment. 

(2)  Costs  of  the  experimental  packages  are  about  the  same  as  Zombie  and 
less  than  the  sub-satellite. 

(3)  The  geometry  is  less  complex  than  for  the  Zombies  and  yields  the  greatest 
possible  separation. 

(4)  The  transmitter  can  operate  for  up  to  six  months  and  thus  may  be  used  with 
a  number  of  receiving  satellites. 

Disadvantages: 

(1)  There  are  only  a  few  P-lls  planned  and  space  has  been  fully  assigned; 
therefore,  it  is  unlikely  that  the  experiment  could  be  flown  unless  another 
experiment  fails  to  materialize  or  is  cancelled. 

(2)  Flights  on  P-11  must  pay  a  weight-proportional  share  of  the  cost  of  the  P-11. 

(3)  RFI  might  be  a  severe  problem  aboard  the  P-11. 
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3.1. 3  Sub-Satellite  to  Agena 


Advantages: 

(1)  The  geometry  is  simpler  than  for  the  other  concepts  and  is  very  precisely 
controlled. 

(2)  The  experiment  Is  less  likely  to  malfunction  than  the  other  experiments. 
Disadvantages: 

(1)  The  weight  is  greater  because  of  the  batteries. 

(2)  The  cost  of  the  experimental  packages  would  be  slightly  higher  than  for 
the  other  experiments. 

(3)  The  experimental  life  of  the  ejected  capsule  would  be  roughly  one  week. 

While  each  configuration  presents  interesting  possibilities,  certain  of  the  practical 
difficulties  or  simplicities  must  dominate  in  the  selection  of  the  optimum  course  of 
action: 

(1)  The  "zombie  to  zombie"  cc  ->.cept  must  be  discarded  for  the  reasons  given  in 
disadvantages  (3)  and  (4)  above.  While  RFI  might  be  eliminated,  the  possi¬ 
bility  of  getting  the  two  packages  into  orbit  within  several  weeks  of  each 
other  is  too  small  to  be  seriously  considered  for  this  type  of  operation. 

(2)  The  "P-11  to  Agena"  concept  must  at  the  present  time  be  discarded  because 
of  disadvantages  (1),  (2),  and  (3)  listed  above.  Again,  while  RFI  might  be 
eliminated,  the  cost  of  the  P-ll  could  be  as  high  as  $250,  000  and,  at 
present,  all  P-ll's  scheduled  to  be  flown  are  completely  filled. 

(3)  The  "sub-satellite  to  Agena"  concept  appears  to  be  the  most  likely 
candidate.  However,  although  this  approach  has  two  very  strong  advan¬ 
tages,  as  listed  above,  the  precise  control  of  experimental  geometry 
would  be  of  particular  importance  in  the  reduction  of  data. 

Therefore  the  "sub-satellite  to  Agena"  concept  is  recommended  for  immediate 
implementation. 
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3.5  EXPERIMENTAL  DESIGN 


The  effects  of  the  ionosphere  upon  the  transmission  between  satellites  will  depend 
upon  the  geometric  relationship  between  the  two  satellites  and  the  ionosphere.  Cal¬ 
culation  of  this  relationship  for  a  large  number  of  points  between  the  start  and  the 
end  of  the  experiment  would  be  expensive  and  wasteful;  therefore,  an  approximate 
technique  was  applied  which  yields  a  rough  solution  above  a  given  altitude.  This  tech¬ 
nique  involves  finding  the  time  from  separation  required  for  the  line  of  sight  (from 
the  mother  vehicle  at  apogee  to  the  sub-satellite)  to  become  tangent  to  an  ionospheric 
layer.  The  mathematical  equations  used  to  calculate  these  data  are  derived  in 
Appendix  A,  which  includes  the  symbols  employed  and  the  Fortran  program  used. 

(It  should  be  noted  that  the  solutions  obtained  are  exact;  the  use  to  which  these  solu¬ 
tions  will  be  put  in  the  following  yields  the  rough  solution  referred  to  above.)  Figure 
3-1,  derived  from  Fig.  A-3  of  Appendix  A,  shows  the  time  required  for  occultation 
of  the  inter-vehicular  line  of  sight  by  the  160-nm  layer  of  the  ionosphere  when  the 
mother  vehicle  is  at  apogee.  This  graph  is  almost  correct  for  the  case  in  which  the 
sub-satellite  is  at  apogee,  because  of  the  small  difference  in  the  orbital  parameters. 

If  it  is  desired  to  find  the  time  of  occultation  for  altitudes  of  the  mother  vehicle  other 
than  apogee,  this  figure  may  be  used  as  follows  to  obtain  an  approximate  answer. 

Given  the  initial  velocity  increment  parallel  to  the  flight  path  (for  an  increment  given 
at  perigee),  enter  the  graph  at  the  time  of  interest  and  read  the  altitude  from  the  apogee 
scale.  This  is  the  rough  altitude  of  the  vehicle  nearest  apogee  (when  both  vehicles  are 
on  the  same  side  of  apogee)  at  the  time  of  occultation  by  the  160-nm  layer.  The  reverse 
procedure  can  be  used  to  find  the  time  of  first  occultation  when  both  vehicles  are  on 
the  same  side  of  apogee  and  the  altitude  of  the  vehicle  nearest  apogee  is  given.  No 
estimate  of  the  error  involved  has  been  made,  but  the  percentage  error  obviously  is 
small  when  both  trajectories  have  a  low  apogee. 

It  is  also  desirable  to  be  able  to  compute  the  inter-satellite  distance  without  resorting 
to  a  computer,  but  it  would  appear  that  this  is  quite  difficult  to  do  since  the  altitude 
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Fig.  3-1  Time  to  Occultation  by  160-nm  Layer  for  Flight  Path  Velocity  Increment  Given 


varies  from  115  nm  to  as  high  as  380  nm.  It  is  shown  in  Appendix  B  however,  that 
the  distance  between  satellite  is  given  within  8  percent  by 


d 


-  2a^a2cos  0^ 


(3.1) 


where 

d  =  distance  between  satellites 
a^  =  semi-major  axis  length  of  first  orbit 
a2  =  semi -major  axis  length  of  second  orbit 
0^  =  earth-centf.i  ed  angle  between  vehicles 

and  the  eccentricities  of  both  orbits  are  under  0. 035;  this  meanB  that  for  a  perigee  of 
115  nm,  the  apogee  is  less  than  245  nm.  Such  an  apogee  is  quite  probable  for  Agena 
flights.  Equation  3.1  may  be  rewritten  as 


d  =  ay  2  -  2  cos  2*  {—  y-)t  (3.1A) 

2  1 


where 

a  =  semi-major  axis  of  the  mother  satellite 
=  period  of  sub-satellite 
T2  =  period  of  mother  satellite 
t  =  time  from  the  first  apogee  following  ejection 

The  reason  for  taking  t ,  above,  from  the  time  of  first  apogee  is  the  fact  that  forward 

ejection  actually  increases  the  orbital  period;  therefore,  the  ejected  package  moves 

ahead  of  and  rises  above  the  mother  vehicle  until  apogee,  at  which  time  it  is  above 

the  mother,  but  the  mother  has  caught  up  and  is  passing  underneath.  For  the  range 

-3 

of  altitudes  under  consideration  here,  the  difference  in  orbital  period  is  1, 5  x  10  hr 
or  5. 4  sec  for  an  initial  separation  velocity  of  7. 5  fps.  At  perigee  at  the  end  of  the 
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first  orbit  following  ejection,  the  interval  between  vehicles  will  be  2.7  sec  or  11.5  nm, 
with  the  mother  now  preceding  the  sub-satellite.  Given  the  7.5  fps  velocity  increment 
at  separation,  the  altitude  difference  at  apogee  will  be  4.7  nm  for  all  altitudes  in  the 
range  under  consideration. 

The  transmitter  structure  will  consist  of  the  transmitter  shell  and  internal  bracketery, 
intermediate  rupport  structure,  and  four  primary  support  fittings.  The  shell  will  be 
attached  to  the  Agena  primary  structure  by  means  of  the  four  primary  fittings  and 
pyrotechnic  pin-puller  assemblies.  The  intermediate  support  structure  is  to  be 
attached  to  the  equipment  rack  with  bolts.  The  release  mechanism  and  impulse  spring 
assembly  will  be  integral  with  the  inte*rmediate  support  structure  and  the  spring  as¬ 
sembly  is  to  be  actuated  with  a  pin  puller. 

Ejection  and  spin  will  be  accomplished  by  pivoting  the  package  about  one  end  and 
releasing  after  90  deg  of  rotation.  A  spring  system  will  produce  an  Initial  tangential 
release  velocity  of  7. 5  fps  and  a  spin  rate  of  1.4  rev/sec.  A  diagram  of  the  ejection- 
spin  up  operation  is  given  in  Fig,  3-2.  The  antennas  are  to  be  erected  at  a  rate  of 
.045  fps,  assuring  an  ample  safety  factor  for  theoretical  loads.  However,  it  is 
planned  to  test  the  antenna  to  determine  its  maximum  moment  of  inertia.  After 
release,  the  package  will  be  spinning  about  the  axis  of  maximum  mass  moment  of 
inertia,  and  due  to  antenna  erection,  the  mass  moment  of  inertia  will  increase  by  a 
factor  of  20.3  in  the  plane  of  rotation.  The  terminal  spin  rate  will  be  .  06  rev/sec. 
These  parameters  are  deemed  satisfactory  for  a  stable  configuration. 

A  preliminary  thermodynamic  analysis  of  the  ejected  package  was  performed  on  the 
basis  of  preliminary  design  data  and  is  detailed  in  appendix  C.  The  mean  vehicle 
temperature  was  calulated  for  a  variety  of  absorptivity-to-emissivity  ratios,  for  two 
spin  orientations,  and  for  several  launch  times.  The  spin  orientation  chosen  <n  the  pre¬ 
liminary  design  (spin  axis  perpendicular  to  earth's  surface  at  perigee)  proves  to  be 
close  to  optimum  for  a  wide  range  of  launch  times.  It  will  be  necessary  to  select 
a  ratio  of  a/e  (absorptivity  to  emissivity)  that  is  suitable  for  all  launch  times; 
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Fig.  3-2  Satellite -to -Satellite  Propagation  Experiment 


but,  for  the  temperature  range  desired  for  the  equipment,  the  mean  temperature 
would  lie  within  a  22  F  range.  A  sample  computer  calculation  showed  that  internal 
major  component  temperatures  will  probably  vary  about  12°F  above  and  below  the 
mean  temperature.  Thus,  the  internal  temperature  may  be  controlled  within 
moderate  bounds  by  the  simple  expedient  of  varying  the  er/e  ratio. 
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Section  4 
SUMMARY 


The  effects  of  the  ionosphere  upon  electromagnetic  radio-frequency  radiation  have  been 
discussed.  Numerical  examples  of  attenuation,  Faraday  effect,  pulse  dispersion,  and 
dispersive  Doppler  have  been  given.  It  was  shown  that  attenuation  would  become  signi¬ 
ficant  at  long  ranges  (1 , 000  nm  or  so) ,  Faraday  effects  would  be  detectable  within  a  few 
minutes,  pulse  shape  deformation  can  be  quite  drastic  at  ranges  of  40  nm  or  less,  and  the 
dispersive  Doppler  effect  is  so  small  as  to  be  virtually  unmeasurable. 

Two  basically  different  approaches  to  ionospheric  studies  were  presented,  the  sounding 
technique  and  a  phenomenological  approach,  and  it  was  indicated  that  the  latter  was 
more  in  keeping  with  satellite  operations.  Several  classes  of  experimental  geometries 
were  discussed,  as  were  the  experiments  possible  with  each  and  some  of  the  more 
serious  attendant  problems. 

Details  of  the  Government-furnished  equipment  were  given,  and  the  three  possible 
modes  of  operation  within  the  limits  of  Program  162  were  described  in  detail.  After 
a  consideration  of  all  factors,  it  was  decided  that  the  experiment  could  be  performed 
only  as  originally  envisioned,  with  ejection  of  a  small,  transmitting  satellite  from  the 
Agena.  Details  of  the  orbital  geometry  were  calculated,  and  an  acceptable  method  of 
achieving  ejection  and  spin  was  designed.  Thermal  calculations  were  performed  to 
assure  that  internal  temperatures  could  be  maintained  so  as  not  to  exceed  equipment 
design  limits. 
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Appendix  A 

LINE-OF-SIGHT  OCCULTATION  BY  ATMOSPHERIC  LAYER 

In  conjunction  with  the  proposed  experiment  investigating  radio  propagation  between 
earth  satellites,  it  was  necessary  to  determine  the  time  required  to  achieve  occulta- 
tion  by  an  ionospheric  layer  of  the  line  of  sight  between  a  mother  and  an  offspring 
satellite.  It  seemed  advisable  to  consider  the  problem  in  its  most  general  form  and 
to  specialize  to  the  particular  problem  only  in  the  last  steps.  Accordingly,  the  rela¬ 
tionship  was  found  between  the  original  orbital  parameters,  an  arbitrary  velocity 
increment  given  to  the  satellite,  and  the  resultant  orbital  parameters  -  the  only  re¬ 
strictive  assumptions  being  that  the  velocity  increment  is  taken  as  impulsive  and  that 
the  earth's  gravity  field  is  taken  to  be  that  resulting  from  a  point  mass.  From  the 
geometry  of  the  system  of  satellites  and  earth,  a  family  of  configurations  was  found 
at  which  occultation  can  occur  -  assuming  the  major  axes  of  the  mother  and  offspring 
orbits  are  coincident  and  that  the  mother  iB  at  apogee  when  occultation  occurB. 

Finally,  a  computer  program  was  written  and  run  to  find  the  time  to  occultation  for 
various  mother  satellite  orbits  using  the  results  from  above  and  assuming  the  entire 
(small)  velocity  increment  to  be  along  the  flight  path.  The  program  computes  the 
time  from  incrementation  until  the  satellites  have  relative  positions  such  that  occul¬ 
tation  will  occur  the  next  time  the  mother  satellite  passes  its  apogee  -  not  necessarily 
the  time  from  incrementation  to  the  first  occultation.  The  parameters  studied  are: 
apogee  and  perigee  height  of  the  mother  satellite;  incrementation  velocity,  and  position 
in  orbit  at  which  the  velocity  is  incremented.  The  computer  output  includes:  mother 
satellite  period;  offspring  period;  mother  satellite  orbit  semi-major  axis,  eccentricity, 
and  apogee  height;  and  the  "Time  to  Occultation." 

LEST  OF  VARIABLES* 

a  semi-major  axis  of  ellipse 

e  eccentricity  of  ellipse 


♦See  Fig.  A-l 
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Fig.  A-l  Orbit  Variables 
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acceleration  due  to  gravity  at  earth's  surface 
special  grouping  of  variables  -Vp/a(l  -  e2) 
dimensionless  velocity  -  v 

6  S ill  6 

special  grouping  of  variables  -  - 5 

i  +  6  COS  u 
2 

focus  to  orbit  distance  -  a(l  -e  )/l  +  e  cos  6 
mean  radius  of  earth 


special  grouping  of  variables 


ll  +  e2  +  2e 
V  1  -  e5 


2e  cos  9 


period  of  satellite  orbit 

velocity  increment  imparted  to  satellite 

velocity  of  satellite  ,  / - ^ J 1  +  e2  +  2e  cos  6 

\l  ,1  2 .  \  o  o  o 

V  a(l  -  e  ) 

angle  between  velocity  increment  and  original  velocity 
angle  between  velocity-velocity  increment  plane  and  orbital  plane 
angle  between  original  and  final  orbital  planes 
angular  distance  from  perigee 

2 

gravitational  constant  of  earth  -  gR 

angle  between  velocity  vector  and  normal  to  radial  line  from  occupied  focus  to 
satellite. 


SUBSCRIPTS 
o  original  orbit 

1  offspring  orbit  or  1-axis  component 

2  2 -axis  component 

3  3 -axis  component 

L  value  of  variable  when  velocity  increment  is  given 

Define  the  1,  2,  3  coordinate  system  as  follows:  the  1  axis  lies  along  the  velocity 
vector  in  the  original  orbit;  the  2  axis  lies  normal  to  1  in  the  plane  of  the  original 
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orbit  with  the  positive  sense  in  the  direction  of  increasing  orbital  altitude;  the  3  axis 
lies  perpendicular  to  the  plane  of  the  orbit  with  positive  sense  such  that  1,  2,  and  3 
form  a  right-handed  coordinate  system. 

ip  is  the  angle  between  the  radial  line  and  the  2  axis-positive  when  the  velocity  has  a 
positive  radial  component. 

Choose  the  angles  a  and  0  such  that  an  arbitrary 
line,  PQ  ,  is  defined  by  them  in  the  following  way: 
a  is  the  angle  between  the  1  axis  and  PQ  - 
0  deg  <  a  <  180  deg  -  and  0  is  the  angle  between 
the  1-PQ  and  the  1-2  planes  - 
-  180  deg  <  a  <  180  deg 


If  a  velocity  increment  v  is  given  to  the  satellite, 
denoting  the  values  of  r  ,  6  ,  V  ,  and  ip  at  this  time  by  the  subscript  L  : 


and 


VL1  +  Vj^  =  VL  +  v  cos  a 

VL2  +  v2  ~  v  8in  a  cos  P 

^L3  +  v3  =  v  sin  a  sin  P 


VL  +  y2  + 


2V^  v  cos  a 


(A.  1) 


The  new  orbital  plane  lies  at  an  angle  y  to  the  original  plane  determined  by 
tan  y  =  V^/V^p  -  where  V^p  is  the  velocity  in  the  1,  2  plane  perpendicular  to 
the  radius,  thus: 


tan  y  - 


_ v  sin  a  sin  0 _ 

cos  if> l  +  v  cos  a  cos  ip^  -  v  sin  a  cos  0  sin  ^ 


(A.  2) 
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The  new  velocity  vector  makes  an  angle  ^  =  arcsin  V^/Vj  with  the  perpendicular 
to  the  radius,  thus 


VT  sin  p.  +  v  cos  a  sin  +  v  sin  a  cos  p  cos 
sin  2  2  779  “* 

(VL  +  V  +  2VLV  008  a* 


(A.  3) 


For  an  orbit  with  semi-major  axis  aQ ,  eccentricity  eQ ,  and  its  perigee  located  at 
eQ  =  o ,  if  n  =  gR2 


tan  ± 
ro 


e  Bin  0 
o  o 


1  +  e  cos  9 
o  o 


(Ref.  4) 


(A.  4) 


and 


V2  « 


a  (1  -  e  ) 
o'  o  ' 


1  +  eQ  +  2eQ  cos  0n  \  (by  implication)  (Ref.  4)  (A.  5) 


Thus ,  since  V 


circ 


/mTT:  (Ref.  4) 


max 


v2  i1 .+  «r 

circ  (1  -  e2) 


1  t  e  y2 

1  -  e  circ 


min 


=  V 


n  i  _ 


circ  2. 

(1  -  e  ) 


1  +  e  circ 


max 


Furthermore: 


“ot1  -  e.2) 

1  +  e  cos  6 
o  o 
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(Ref.  4) 


(A.  6) 


When  r  .  V .  and  ^  are  known,  the  orbital  parameters  a .  e ,  and  0  may  be  found 
from: 


a  = 


2  -  rV  /n 


from 


v2-M 


(Ref.  4)  (A.  7) 


cob*  ^  (Ref.  3  or  4)  (A.  8) 

(by  implication) 


and 


cot  0  =  cot  ^  |1  2  2 

rV  cos^  <p 


) 


(Ref.  2)  (A. 9) 


Thus,  if  the  velocity  increment  v  is  given  at  0  =  0^,  using  relations  (A.4andA.5) 
with  Eqs.  (A.  1,  2,  and  3)  gives: 


=  K2  (l  +  e2  +  2e  cos  9 .  ) 
o  \  o  o  L' 


+  2e  cos  9.  /  +  v  +  2vK  cos 
o  L'  o 


a  ]/l 


1  +  e  +  2e  cos  9r 
o  o  L 


(A.  10) 


where 


(A.  11) 


while 


tan  y 


_ v  sin  a  sin  B  _ _ 

v  cos  a(I  +  eQ  cos  6^)  eT v  sin  a  cos  p  sin  9 L 


K  (  1  +  e  \  +  —  =====  ==«=  -  .  . = .  i 

0  °  y  1  +  e2  +  2e  cos  9.  y  1  +  e2  +  2e  cos  9, 

7  o  o  L  »  oo  1 


(A.  12) 
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(A.  16) 


2 

v  a. 


1  - 


2  cos  av 


f  \  1 1  +  eo  +  2 

V  <*\ 


2e  cos  0, 
o  L 


Defining 


PT  =  v 


K  riKZ 

V  n  eQ] 

/  eo  sin  flL  \ 

V  +  eo  cos  °L/ 


rp  =  perigee  radius 


S,  s. 


1  +  e  +  2e^  cos  0r 
o  o  L 


1  -  e‘ 


Relations  (A.  14)  and  (A.  15)  become 


\  M  /  1  +  eo  COS  0L  l  \ 


pl  +  2Plsl  008 


“>] 


and 


2  2 
r  V,  cos  di 
o  l 


-  =  (1  +  eQ  cos  0L) 


1  +  2Pt 


cos  a  -  qL  sin  a  cos  0 


+  p: 


2  2 
1  +  q*  -  (qL  cos  a  +  sin  a  cos  0) 


s: 
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(A.  17) 

(A.  18) 

(A.  19) 

(A.  20) 

(A.  21) 


Thus 


(PL  + 


2PlSl  cob 


4 


(A.  22) 


The  expressions  for  e^  and  0^  are  not  as  near  and  simple.  However,  by  calculating 
relations  (A.  20)  and  (A.  21)  by  use  of  Eqs.  (A.  17,  18,  and  19)  and  computing  cot  ^ 
by  use  of  Eq.  (A.  4),  the  values  of  e1  and  0^  may  be  readily  determined. 

Assuming  that  eQ  and  Pj  are  small  leads  to 

qT  «  e  sin  0. 

L  o  L 

S.  »  1  +  e  cos  0T 
L  o  L 

Thus 

(2 

and 

..2  2  . 

r  V*  cos  <p.  , 

- -  «  (1  +  eQ  cos  0L)  11  +  2PL(cos  a  -  eQ  sin  a  cos  0  sin  0^)(1  -  eQ  cos  0 L) 

r 

Thus 

a, /a  «  1  +  2P.  cos  0(1  +  e„  cos  0T  ) 

1  O  Li  O  Li 


4) 


(1 


e  cos 


V 


[>- 


2PT  cos  a(l  +  e  cos  0. 


) 


This  equation  was  derived  by  using  an  approximation  of  actual  e^,  assuming  eQ  and 

p  to  be  small.  The  value  of  0,  would  probably  be  much  different  from  that  predicted 
o 

by  a  small-value  approximation  for  eQ  and  Pq  and  for  that  reason  will  not  be 
approximated  here. 
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The  period  of  the  satellite  is  T  =  2ir 


Thus 


Tl/To  =  >/  <al/ao)3  ~  (l  -PL'  2PLSL  COS  °) 


-3/2 


1  +  Jv  \  |  —  cos  a 

o  \  p 


cos  a 


M 


— 2 - \r3/2 

1  *  %  +  co»  *L  j 

>  -  %  /J 


1  +  e  +  2e  cos  0. 
_ o _ o _ L 

1  -e2 


(A.  23) 


Thus,  if  all  other  variables  are  fixed,  the  greatest  change  in  period  will  arise  from 
giving  the  velocity  increment  at  0Q  =  0  (perigee),  the  least  from  giving  the  increment 
at  6  -  7r  (apogee). 

It  is  interesting  to  note  the  small  effect  on  the  orbital  parameters  due  to  the  component 
of  the  incremental  velocity  which  does  not  lie  parallel  to  the  flight  path.  An  intuitive 
explanation  of  this  is  that  -  so  far  as  the  component  normal  to  the  orbit  is  concerned  — 
the  gravity  force  vector  has  a  component  which  points 
in  such  a  direction  that  the  satellite  will  be  accelerated 
back  toward  the  original  orbital  plane.  The  force  due 
to  this  component  is  (nearly)  proportional  to  the  distance 
from  the  original  orbital  plane  so  simple  harmonic  motion 
occurs  about  the  original  orbital  position  and  is  super¬ 
imposed  on  the  old  orbital  motion.  The  oscillation  has  a 
period  nearly  equal  to  the  orbital  period  so  that  the 
orbital  plane  is  rotated  but  other  effects  are  quite  small. 

For  a  radial  component  -  the  difference  p/r  -  v  /r  (the  downward  force)  decreases 
as  the  satellite  drops  (assuming  v  =  constant )  and  increases  as  the  satellite  rises 
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causing  an  acceleration  back  toward  the 
original  orbit.  The  period  of  the  oscillation 
is  nearly  equal  to  the  orbital  period.  Thus 
a  radial  increment  causes  a  change  in 
eccentricity  with  a  small  change  in  period. 


Defining  pQ 
radii, 


I©  *•  j 

—  ;  p.  =  —  the  normalized 
rf  1  rf 


p  -  sec  (b 
Ho  ^o 

PL  =  sec 

Thus 


<pc  =  arc  sec  pQ  +  arc  sect  =  arc  cos  l/pQ  +  arc  cos  l/p1 


U  <p  is  the  angle  from  the  line  of  apsides  of  the  original  orbit,  9  is  measured  from 
c 

the  original  orbit  perigee,  and  is  the  angle  between  the  lines  of  apsides  of  the 
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two  orbits:  (0  -  ^1)crlt  "  *  +  0c  -  0  and 

,  ai(‘  -  ei) _ i 

pl  R  +  hj  1  -  cos  (tpc  -  Jjj) 

=  Px  1/1  -  ej  cos  (<(>c  -  0X) 

when 


P 


1 


R  +  hf 


(A.  25) 


now 


cos  <j>  =  cos 


But  from  trigonometry: 


arc  cos  l/n  +  arc  cos 


1  -  e  cos  ( <(>  -  0,) 


] 


cos  (a  +  b)  =  cos  a  cos  b  -  sin  a  sin  b 


Thus 


cos  <p  = 


1  -  e,  cos  (<p 


-  v  Vpq  -  1  VP1  -  (X  -  el  COS  K  "  °ll) 


C  p 

“r 


P0P!  008  4>c  +  ej  COB  (4>c  "  V  ”  1 


il  - 1  \[pi  -  (i  -  ei  co>  |+c  -  ®i|)z 


(A.  26) 
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Since  solution  of  this  for  general  0 ^  would  be  very  difficult  and  since  in  the  problem 
under  consideration  is  almost  zero,  it  will  be  assumed  to  be  so,  then 


(poPl  +  ei)cos  *c  ‘  1  =  -Jp o  ‘  1  /Pi  ‘  (J  ‘  ei  008  ^c)2  (A,2?) 

from  which 

[(»„P1  +  *lf  *  el(po  •  ')]co‘2*c  -  2  [(<’„P.  *  e])  +  "l(<>0  ■  I)]“*c 

*  [*  -  Ul  -  i)(’p?  -  *)]  - 0 


thus 


cos  <f>  = 


\  1  2  A 

(poPl 

+  el 

)  +  ellpo  "  0 

\2  2/2  A 

(poPl 

+el 

+  el  (»o  -  l) 

kp,  *  »)  1  -  (po~  ‘)  (p2r') 


(A.  28) 


The  plus  sign  corresponds  to  the  case  when  both  radial  lines  are  on  the  same  side  of  the 
point  of  tangency  so  the  minus  sign  gives  the  answer  desired. 


The  time  equation  is 


t'to  = 


=  g  jarcsin 


2  sin  6 
1  +  e  cos  0 


)-  e^l  -  e2  IT 


sin  0 


e  cos 


20 


10 
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Thus,  since  0  =  *  +  <t>c  ~  and  T 


27T 


arc  sin 


cos  (<*>c  - 

(interpreting  the  arc  sin  as -180°  <  arcsin  <  180”) 


[ — T  ■l,»(*c’91) 

+  ejy1-  en-e,  COS  (<*>  -  0n)| 
1  '^c  1' 


(A.  29) 


where  t  /Tj  is  the  fraction  of  the  orbit  through  which  the  satellite  must  pass  in  going 
from  6  -  ir  -  0  ^  to  0  =  ir  +  <|>c  -  0^ .  The  number  of  hours  after  separation  necessary 
for  the  offspring  to  have  gained  (or  lost)  this  "distance"  relative  to  the  mother  is  given  by 


N 


t 

TjT  (orbits) 


(orbits\ 
hour  / 


(A.  30) 


THE  COMPUTER  PROGRAM 


The  program  utilizes  the  results  of  the  analytical  sections  and  is  based  on  the  following 
important  assumptions: 

1.  The  initial  orbit  is  of  small  eccentricity. 

2.  The  incrementation  occurs  near  perigee  or  near  apogee  —  i.e. ,  not  near 
0  =  90°  or  270’ . 

3.  The  ratio  of  the  incremental  velocity  to  the  orbital  velocity  is  small. 

4 .  The  component  of  the  incremental  velocity  normal  to  the  orbital  flight  path 
has  only  a  negligible  effect. 

5.  The  line  of  apsides  (major  axis)  of  the  orbit  is  not  shifted  by  the  incrementation. 

6.  The  occultation  occurs  with  the  mother  satellite  at  apogee. 

7.  The  occulting  layer  is  below  the  apogee  heights  of  the  satellite  orbits. 

8.  "Time  to  Occultation"  is  the  time  from  incrementation  to  the  state  at  which 
occultation  will  occur  the  next  time  the  mother  satellite  is  at  apogee . 

*The  orbital  geometry  is  given  in  Fig.  A-2  and  the  results  of  the  calculation  in  Fig.  A-3. 
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BETWEEN  THE  SATELLITES  PASSES  THROUGH  THE 
IONOSPHERIC  LAYER.  IT  WAS  ASSUMED  THAT  THE 
MOTHER  SATELLITE  IS  AT  APOGEE  WHEN  OCCULTATION 
OCCURS. 


QN033S/133J  -1N3W3M3NI  AXI3013A  HIVd  IHOHJ  0NO1V 
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To  use  the  computer  program  provide  ALH  (see  the  list  of  Fortran  variables),  PH, 
AHI,  AH(  1 ),  TH(  1 ),  V(l),  THI,  VI,  KA.  KV,  KT. 


Fortran 

A(K) 

AH(K) 

AHI 

ALH 

B 

C 

CE 

D 

E(K) 

EE 

EEE 

G 

H 

J 

K 

KA 


ALL  DIMENSIONS  IN  NAUTICAL  MILES 
ALL  SPEEDS  IN  FEET  PER  SECOND 
ALL  ANGLES  IN  RADIANS. 

Variables 

Semi-major  axis  of  mother  orbit  =  aQ 
Apogee  height  of  mother  orbit  =  h 

a 

Difference  between  successive  mother  orbit  apogee 
Altitude  of  occulting  atmospheric  layer  =  h^ 
Semi-major  axis  of  offspring  orbit  =  a1 
Differential  cyclic  frequency -mother  to  offspring 
Cos  (  E  -  7T  ) 

Eccentricity  of  offspring  orbit  =  e^ 

Eccentricity  of  mother  orbit  = 

Arc  tan  (sin  E  -  w/cos  E  -  ir) 

E  -  7T 

Gravitational  constant  =  n 

Approximate  time  to  occultation 

Index  variable  on  location  in  orbit  at  which  velocity 

Index  variable  on  mother  orbit 

Upper  limit  on  value  of  K  -  thus  number  of  orbits  i 


(nautical  miles) 
(nautical  miles) 
heights 

(nautical  miles) 
(nautical  miles) 
(orbits  per  hour) 


(hours) 

is  incremented 

ip  to  25 
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KT  Upper  limit  on  value  of  J  -  thus  number  of  locations  for  incrementation 

up  to  6 

KV  Upper  limit  on  value  of  L  -  thus  number  of  incremental  velocities  up  to  20 

L  Index  variable  on  velocity  increment 

2 

P  Special  constant  =  Pj  =  a^l  -  ej)/(Re+  h^) 

PH  Perigee  Height  =  h 

12  2  1 
(PDPi  +  ei)  +  Zj 

2  2 

QA  Special  constant  =  +  e^)  +  e^  Z 

R  Special  constant  =  cos 

Ro  Special  constant  =  pQ  =  Rg  +  ha/Re  +  h^ 

S  Sin  E  -  7r 

TB(J,K,  L)  Period  of  offspring  orbit  =  Tj  (hours) 

TH(J)  Angle  where  increment  is  given  =  0j 
THI  Difference  between  successive  values  of  0  ^ 

TM(K)  Period  of  mother  orbit  =  Tq  (hours) 

V(L)  Velocity  increment  along  flight  path  =  v  cos  a  (fps) 

VI  Difference  between  successive  velocity  increments 


X 


Y 


Special  function  =  2 


'i 


1  +  e  +  2e  cos  0. 
oo  1 


Inverse  of  circular  velocity 


Z 


Special  constant  =  z  =  p  -  1 
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this  I  S  A  PROGRAM  FOR  UE  TE  RM  I N  ING  THE  PER  10Q_OF  AN  OFFSPRING _ 

SATELLITE  ANC  THE  TIRE  NECESSARY  FOR  OCCULTATICN  OF  I HE  UNE-OF-SIGMI 
PEI  RE  EN  [Hfc  MOTHER  AMI  OFFSPRING  SATELUIES  TO  OCCUR.  I  HE 
PARAMETERS  STUnien  ARE  THE  ORIGINAL  ORBIT,  THE  AL  II  TUOE  OF  THE 
OCCLUOING  LAYER,  THE  ALONG-FL  IGHT-PATH  VELOCITY  INCREMENT  GIVEN  fo 
TO  THE  OFFSPRING,  ANU  IHCLOCAIION  IN  THE  ORIGINAL  CRM  1  T  AT  WHICH  THE 
INCREMENI  IS  GIVEN.  THIS  NEGLECTS  AiR  DRAG,  OHLAIINtS'S  CF  THE 
FARTH,  ANO  ASSUMES  THAT  THE  VELOCITY  INCREMENT  IS  SMALL  AND  GIVEN 
IMPULSIVELY  ANO  THAT  THE  ORIGINAL  ORBIT  IS  OF  SMALL  ECC E N I R I C  1 1 Y . 

DIMENSION  THI6I , VI 20  I ,AH(2S» ,AI 251 ,E (251 , IMI25 I ,  Tit  <2  5, 6, 20) ,HI25,6 

7 , 20 ) ,N ( 6 ) 

READ  INPUT  TAPE  5 ,  1 0 , ALH, PH, TH I  II , AH  I  I  I , V II) , TH I , All  I  ,  V  I , KA , KV ,KT 
10  FORMAT (4E 15.4 /ME  IS. 4 ,4  1ST 

WR  ITEOUTPUTTAPE  6,  10, ALH, PH, TH I  I  ) ,  AH (  II , V (  I  I , T H I , AH I , V  I , K A , KV , KT 
G«i2.  I  /  ■  (  3960.* 5 2H0.  1--2/6076. I 
20  0095K-I.KA 

AIK)«3440.*TAHIK JfPHI/2. 

_ E(K)»( AHtK)-Phl/l2.-3440.*AH|KT»PHI _ _ _  _ 

R0«(3440.*AH|kF|/(  3440. ♦ALH) 

Z«RO»«2-  I. 

Y-SORTF(A(K)/G) 

TM(K)«2.«3.l4l5927*Y»A(Kl-6.0761/3.6 
30  0090 J* I ,KT 

X-SORTFI (  1.*6<KI«(EIKT+2.»C0SF ITHI Jill  I/I  I  .-E ( K I »»2 I ) « V-2. _ 


NO  0080L* I , KV 

TB(K,J,L)«TMIK»/m.-V(L  I.S)  _ 

B*  A I K ) / ( I .-V I  L  I  *X  1 

C-U  1.-V(LT»X»»»I.5-1.»/TNIKI  __ 

O-SQRTF IE(K l««2»4.»V(L )»Y»E (K» -C0SFI TH<J>l*4.-(V<L)»Y)*»2t 

P-B» (J .-0«»2J /<  3440.»ALH> _ 

0A»rR0»PF0l»»2FZ«n*»2 
0*  (  RO*  P  ♦  D«  (  1  .  *  /  I  l/OA 
R-Q-SQRTF IQ--2-I 1 .  —  Z • C P-»2- I . ) I/CAI 
S-SQRTFI ( I .-0--2TM 1 .-R»»2) )/(  l.-0«K> 

CE » C  0-K 1 / I  I •  —  0 • R  » 

EE-AIANFIS/CEI 


IF  I  EE  I  SO , 50 , 60 
50  EEE— EE 
GO  TO  70 

60  EEE-3. 1415927-EE 

70  HI  K,  J,L  )"IEEE-»S»DI/I6.?B3  IB5  3»ABSF  (C  >  I 
VIL*H»V(LT«VI 


IF  IV IL ♦  1 1  TWO , 75,80 
75  VIL+H-VIL  1*2. *VI 
80  CONI INUE 
90  TH(J»I)«THI J)«TH 
95  AHIKFl I  « AH ( K  I  HAH  I 

_ 00  1001-  I  ,K  T _ _ _ _ 

100  N ( I  I - 1 

COI90K-  l,KA  . 

WRITE0UTPUTTAPE6, IT 0 , K , AHl K ) ,K , A ( K ) , K , TM I K I , K, E I K I 
T  10  FORMAT!  IH8,UH  Ah  I  1 2 , 2H I -£ I  2. 4 , 5X ,2HA I  12  * 2H I -E 1 2 . 4 , 5X , 3HTMI 12, 2HI »E 
.  712.4,5X,2HE( I2.2HI-E 12.4) 

_ WKJ TEOUTPUTTAP E_6j  I2C  .  I N(  I  )  ,  1*1  ,  K  T ) _ _ 

120  FORMAT  I  1  HO ,  I0X , 6 1  I  IX.3HTHI I  I. IMI II 
WR  ITE0UTPUTTAPE6, I  30 , I THIII , I- I , KT  | 


THIS  IS  A  PROGRAM  FOR  DETERMINING  THE  PERIOD  OF  AN  OFFSPRING 

130  FORMATC  IMO,  I7X, IE I2.R.RX) T 

WRIT  ECU!  PU  T  TAP'  o  ,  I  :.u  ,  (  X. ,  N  (I  T  ,  I  *  I  ,K  T  ) 

I  MO  FORMAT  I  2H0L  i.tVIL  I  ,2X,6(6X,  SHI  Hi  12.  IH,  II  ,  3H,L  I  I  I 

no  isoi  i./.v 

ISO  WR 1TCOUI  PUTT  APE  6, IA0.L,V(LI,1I!»IK,J,LI.J* I  »KT  > _ 

I  AO  FORMAT!  IH  ,  1 2 ,  I X , E 1 2 . R , 2X , E 1 2. R , I RX , E 1 2 . R I ) 

WR  I  T  EOll  f  PU  T  T  APE  6  ,  I  70  ,  I K  ,N  (  IT  ,  I  *  T  ,K  I  » 

TIO  FORMAT (2H0L,5X»RHV(L  T .  2X  ,  6  (  7  X ,  2HH  (  I  2 ,  ITT*  I  I  ,  SM,U  II 
00I80L- I.KV 

IPO  WR I TEOUI PUT  TAPE  6  » I6G »L . V ( L I , (M I K, J, L » , J* 1 ,KT I 

190  CONTINUE _ _ _  _ _ 

"CALL  EXIT 

ENCI  1, 0,0, 0,0.0. 1.0,0, 1,0,0,0,0,01 
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NOTE:  The  equations  which  refer  to  the  sources  may  all  be  derived  from  the  basic 
equations 


%(>  -  °2) 

1  +  e  cos  0 


and 


tan  ^ 


I 

r  d  9 


e  sin  0 
1  +  e  cos  0 


all  of  which  may  be  found  in  Ref.  4. 
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Appendix  B 

APPROXIMATION  TO  INTER-SATELLITE  DISTANCE 


If  one  knows  the  distance  from  each  satellite 
to  the  center  of  the  earth  and  the  angle  between 
the  lines  passing  from  the  earth  to  the  satellites, 
the  distance  between  the  satellites  is  found  by 
the  law  of  cosines  to  be 

d2  =  ri  +  v  2rir2  cos  9d 

using  the  series  for  the  inverse  of  the  time 
equation  as  given  by  Moulton,  et  al. ,  in 
Periodic  Orbits,  Washington,  1920  (pp.  55-59). 
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Thus 


Since 


if 


For 


thus 


2  2  2 
ri  =  ai  +  ai 


-  2et  cos  cjj  t  +  (1  -  cos  2u)jt) 


3  2  2 

-  cos  0)^(1  -  cos  2o^t)  +  cos  Wjt 


iv/  ~i 


i  2 

+  —r  (1  -  COS  2u>  t)  +  .  .  . 


cos  (a  +  b)  =  cos  a  cos  b  -  sin  a  sin  b 


fl(t)  s  (2e2  -  -4  ]  sin 


/  °i\ 

w2t  -  |2ex  -  -4-  Isin  Wjl 


5e^  5eZ 

+  sin  2w^t - sin  2u>jt  + 


et  w  0.035  ,  0(t)  =s  0.10 


cos  O(t)  »  1.0 
sin  O(t)  »  fl(t) 
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Hence 


cos  0d  =  cos  (w2  -  w1)  t  -  O(t)  sin  (o>2  -  w  )  t 


and 


2  2.2 


d"  = 


ax  +  a2  -  28^2  cos  (w2  -  t  +  2a1a2f2(t)  sin  (o>2  -  Wj)  t 

2  (  2  2  2  ) 
+  aj^  <-  2e^  cos  a>1 1  i  e^(  1  -  cos  2w^  t)  +  cos  t  +  . . . .  | 

2  |  2  2  2  ) 
+  a2  2e2  cos  w2t  +  e2(l  -  cos  2u>2t)  +  e2  cos  w2t  + - J 


'ia2  [ 


-  2a, a,  |  cos  (w2  -  )  t  -  0(t)  sin  (u>2  -  c^)  t 


-  ej^  cos  Wjt  -  e2  cos  w2t 


2  2 

ei  e2 

+  e^  cos  u^t  cos  Wgt  +  -y  (1  -  cos  2w1t)  +  -g-  (1  -  cos  2u>2t)  +  .... 


If 


.2  2  2  2 

d  =  ax  +  a2  -  2a1a2  cos  (w2  -  u^)  t  +  t  a^ 


the  error  term  e  would  not  be  expected  to  be  larger  than 


M 


=  +  4{c2  -  ex)  +  §(e*  -  ej)  +  {-  2c  x  +  3c*J  +  {-  2e2  +  3e*J 

2  2] 

3le2  +  el  +  e2 J 


“  el  "  ®2  +  e’ 


»  -  8e, 
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when 


COS  (u>2  -  Idj )  t  =  0 


then 


,2  2  2  2 

d  =  at  +  a2  +  c  &1  a2 


if 


al  =  a2  *  rf2  =  <2  +  e2)a2=  1.72a2 


2  2 

and  Bince  d  «  2a 


then 


d>ctual  *>  1'3Ia 


d  »  1.41a 

approx 


hence,  the  error  is  8  percent  or  less  if  one  says 


d2  =  ai  +  a2  -  2ai  a2  008  9d 


or 


=  yjn*  +  a2  "  ^1  a2  008  <w2  '  "l^  ’  if  ei  - 
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e2  so.  035 


Appendix  C 

PRELIMINARY  THERMODYNAMIC  ANALYSIS 


The  vehicle  was  assumed  to  be  a  box  of  .  04-in.  thick  HM214-T8  magnesium  10  by  14  by 
20  in.  covering  several  internal  equipment  boxes.  Total  weight  was  assumed  to  be  50  lb, 
of  which  5  lb  was  external  structure.  The  equipment  was  assumed  to  have  2.6  w  contin¬ 
uous  power  dissipation  internally.  The  experiment  was  assumed  injected  into  an  orbit 
between  0  ~  0  and  (3  =  45  deg  (60  to  64  percent  sunlight)  and  to  be  spinning  about  its 
stable  axis.  The  orbit  was  assumed  to  vary  from  115  to  225  nm  and  was  nearly  polar. 
Calculations  were  performed  for  the  cases  of  spin  axis  at  parallel  and  spin  axis  perpen¬ 
dicular  to  the  earth's  surface  when  the  satellite  is  perigee. 

An  analytical  approximation  of  the  mean  orbital  temperatures  for  the  entire  vehicle 
may  be  made  utilizing  time-averaged  heat  rates  for  the  external  surfaces  and  using 
different  paint  patterns  with  varying  ratios. 

The  differential  equation  for  satellite  temperature  (assuming  the  spin  axis  parallel  to 
the  earth  surface  at  perigee)  is: 

c~  =  £  AofS  +  £  AoEa  +  £  AeEs  +  P  -  £  AeoT4  (C.  1) 


where 

c  =  thermal  capacity  ~  Btu/‘  F 

T  =  mean  satellite  temperature  ~  °  F 

2 

A  =  external  surface  area  ~  ft 
a  -  absorptivity 
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f  =  fraction  of  orbit  in  sunlight 

2 

S  =  solar  constant  ~  .  123  Btu/sec-ft 

e  =  emissivity 

Eg  =  earth  radiation  in  infrared  spectrum 

E  =  earth  reflection  or  albedo 
a 

P  =  power 

cr  =  Stefan- Boltzman  constant 

9  =  time 

2 

The  maximum  area  that  can  be  illuminated  by  the  sun  is  2.456  ft  and  minimum  is 
1.505  ft2. 

2  2 
The  leading  surface  in  orbit  is  1.945  ft  and  has  average  earth  radiation  of  .007  Btu/ft  - 

2  2 

sec  and  average  albedo  of  .  00515  Btu/ft  -sec  in  noon  orbit  and  .  00364  Btu/ft  -sec  in 

j3  =  -45  deg  orbit.  These  heat  rates  are  symmetrical  for  the  trailing  edge.  The  circum- 

2  2 
ference  has  4.73  ft  of  surface  area  and  average  earth  radiation  of  .0068  Btu/ft  -sec 

2  2 
and  albedo  of  .  0049  Btu/ft  -sec  in  a  noon  orbit  with  a  value  of  .  00347  Btu/ft  -sec  in 

P  =  -45  deg  orbit. 

In  order  to  simplify  the  solution,  it  will  be  assumed  that  the  outer  paint  surface 

a/e  ratio  is  varied  for  p*  =  0  and  p  =  45  deg.  In  p  =  0  orbits,  60  percent  of  orbit 

is  in  sunlight  and  64  percent  is  in  sunlight  for  p  =  -45  deg.  Substituting  values  in  the 

H  T 

equation,  remembering  that  c  -j-^  =  0  in  steady  state,  and  making  suitable  modification 
to  Eq.  (C.  1)  gives: 

P  =  0  deg  (MINIMUM  ENERGY  ORBIT,  spin  axis  parallel  to  earth  at  perigee) 

<rT4  =  .0179  a/e  +  .00689  (C2> 


*p  =  angle  between  earth-sun  line  and  satellite  orbital  plane. 
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(C.3) 


p  =  ±  45  deg  (spin  axis  parallel  to  earth  at  perigee) 
o-T4  =  .  02597  a/e  +  .  00689 

p  =  0  to  p  ='  *  90  (spin  axis  passes  through  earth  center  at  perigee) 

4  _  ( 1, 945  cos  P  +  1.  505  sin  P  .  123f  +  .  005  cos 
aT  ~  \  8.62 

Equations  (C.  1)  and  (C.2)  are  plotted  in  Fig.  C-l  and  various  a/e  values  for  Eq.  (C.3) 
are  plotted  in  Fig.  C-2. 

In  order  to  get  typical  time-temperature  histories  of  the  satellite,  a  small  computer 
mathematical  model  was  set  up  utilizing  4  nodes  or  equipment  parts.  Node  1  cor¬ 
responded  to  the  leading  edge  at  perigee  and  Node  3  to  the  trailing  edge  at  perigee. 

Node  2  represented  the  other  external  surfaces.  Node  4  was  the  lumped  interior  equip¬ 
ment,  including  the  transmitter. 

Orbital  heating  rates  were  entered  as  time-dependent  functions  into  external  surfaces 
to  account  for  solar  illumination,  earth  reflection,  and  earth  emission  in  the  infrared 
spectrum.  The  internal  power  was  delivered  to  the  internal  equipment  node. 

The  vehicle  interior  was  assumed  to  be  one  homogeneous  structure  of  similar  dimensions. 
External  structure,  internal  surface,  and  internal  box  were  painted  with  black  Kemacryl  to 
maximize  internal  radiative  coupling.  This  will  permit  the  2.6  w  to  be  transmitted  to 
vehicle  skin  and  to  be  rejected  to  space.  This  analysis  assumed  further  that  there  was 
no  convective  or  conductive  heat  transfer  from  vehicle  interior  to  the  surface.  An 
a/e  ratio  of  1. 1:1  was  used. 

The  most  likely  orbit  is  expected  to  be  that  where  p  =  -45  deg  with  the  solar  vector  in¬ 
clined  38. 5  deg  to  the  spin  axis.  This  was  the  case  selected  for  analysis  on  the  digital 
computer.  The  mathematical  model  then  produced  a  time- temperature  history  of  the 
equipment  node  and  skin  nodes  showing  the  effects  of  entering  and  leaving  earth's  shadow 
and  other  orbital  variations  experienced  by  the  vehicle. 


|  a/e  +  .00689  (C.4) 
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as  a  Function  of  Other  P  's 


C-5 


CONCLUSIONS 


If  the  following  conditions  are  furnished  or  established,  the  vehicle  internal  components 
can  probably  be  maintained  within  temperature  limits.  First,  excellent  radiant  coupling 
of  interior  components  with  external  surfaces  must  be  maintained.  This  is  easily  accom¬ 
plished  by  selection  of  a  high-emissivity  paint.  Also,  with  the  expected  complex  internal 
structure,  conductive  heat  transfer  must  be  maximized  so  all  component  mountings  should 
be  carefully  covered  with  a  conducting  silicone-type  grease. 

The  internal  power  dissipation  and  sequence  of  events  must  be  known  very  accurately. 

A  mandatory  requirement  is  precise  definition  of  launch  window,  ejection  time,  dynamic 
history  and  final  sun/satellite  geometrical  relationships. 

It  is  possible  to  prescribe  external  finishes  on  the  experiment  which  will  maintain  tem¬ 
peratures  within  the  specified  limits  for  orbital  lifetimes  on  the  order  of  2  weeks.  How¬ 
ever,  due  to  many  uncertainties  and  changes  in  energy  fluxes,  a  passive  system  will  not 
successfully  operate  over  3-month  or  longer  lifetimes. 
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DISTRIBUTION 


Code 
AF  5 

AF  18 

AF  32 

AF  33 

AF  43 
AF  124 

AF  139 
AF  314 
AF  318 

Ar  5 


Ar  9 


Ar  50 


No.  of  Copies 


AFMTC  (AFMTC  Tech  Library-MU-135)  1 

Patrick  AFB,  Fla. 

AUL  1 

Maxwell  AFB,  Ala. 

OAR  (BROS,  Col.  John  R.  Fowler)  1 

Tempo  D 

4th  and  Independence  Ave,  Wash  25,  D.  C. 

AF06R,  OAR(SRYP)  1 

Tempo  D 

4th  and  Independence  Ave,  Wash  25,  D.  C. 

ASD  (ASA  PRD  -  Dist)  1 

Wright- Patterson  AFB,  Ohio 

RADC  (RAYLD)  1 

Griffiss  AFB,  New  York 
Attn:  Documents  Library 

AF  Missile  Development  Center  (MDGRT)  1 

Holloman  AFB,  New  Mexico 

Hq.  OAR  (RROSP,  MaJ.  Richard  W.  Nelson)  1 

Washington  25,  D.  C. 

ARL  (ARA-2)  1 

Library  AFL  2292,  Building  450 
Wright- Patterson  AFB,  Ohio 

Commanding  General  1 


USASRDL 

Ft.  Monmouth,  N.  J. 

Attn:  Tech  Doc.  Ctr. 

SIGRA/SL-ADT 

Department  of  the  Army  1 

Office  of  the  Chief  Signal  Officer 
Washington  25,  D.  C. 

Attn:  SIGRD-4a-2 

Commanding  Officer  1 

Attn:  ORDTL-012 

Diamond  Ordnance  Fuze  Laboratories 
Washington  25,  D.  C. 


D-l 


Code 

Organization 

No.  of  Copies 

Ar  67 

Redstone  Scientific  Information  Center 

U.  S.  Army  Missile  Command 

Redstone  Arsenal,  Alabama 

1 

O  31 

Office  of  Scientific  Intelligence 

Central  Intelligence  Agency 

2430  E  Street,  N.  W. 

Washington  25,  D.  C. 

1 

Q  2 

ASTIA  (TIPAA) 

Arlington  Hall  Station 

Arlington  12,  Cirginia 

20 

G  109 

Director 

Langley  Research  Center 

National  Aeronautics  and  Space  Administration 
Langley  Field,  Virginia 

1 

N  9 

Chief,  Bureau  of  Naval  Weapons 

Department  of  the  Navy 

Washington  25,  D.  C. 

Attn:  D LI-31 

2 

N  29 

Director  (Code  2027) 

U.S.  Naval  Research  Laboratory 

Washington  25,  D.  C. 

2 

I  292 

Director,  USAF  Project  RAND 

The  Rand  Corporation 

1700  Main  Street 

Santa  Monica,  California 

THRU:  AF  Liaison  Office 

1 

M  6 

AFCRL,  OAR  (CRXRA  -  Stop  39) 

L.  G.  Hanscom  Field 

Bedford,  Mass,  (mail  separately) 

20 

AF  253 

Technical  Information  Office 

European  Office,  Aerospace  Research 

Shell  Building,  47  Cantersteen 

Brussels,  Belgium 

1 

Ar  107 

U.S. Army  Aviation  Human  Research  Unit 

U.S.  Continental  Army  Command 

P.  O.  Box  428,  Fort  Rucker,  Alabama 

Attn:  Maj.  Arne  H.  Eliasson 

1 

D-2 

No.  of  Copies 


Code  Organization 


G  8  Library  2 

Boulder  Laboratories 
National  Bureau  of  Standards 
Boulder,  Colorado 

M  63  Institute  of  the  Aerospace  Sciences,  Inc.  1 

2  East  64th  Street 
New  York  21,  New  York 
Attn:  Librarian 

M  84  AFCRL,  OAR  (CRXR,  J.  R.  Marple)  1 

L.  G.  Hanscom  Field 
Bedford,  Mass. 

N  73  Office  of  Naval  Research  10 

Brance  Office,  London 
Navy  100,  Box  39 
F.  P.  O. ,  New  York,  N.  Y. 

U  32  Massachusetts  Institute  of  Technology  1 


Research  Laboratory 
Building  26,  Room  327 
Cambridge  39,  Massachusetts 
Attn:  John  H.  Hewitt 


IT  431  Alderman  Library  1 

University  of  Virginia 
Charlottesville,  Virginia 

G  9  Defence  Research  Member  1 

Canadian  Joint  Staff 
2450  Massachusetts  Avenue,  N.W. 

Washington  8,  D.C. 

AF  353  ESD  (ESRDC,  Mr.  Jonas  Hallgrimson)  1 

L.  G.  Hanscom  Field 
Bedford,  Mass. 

Lockheed  Missiles  &  Space  Co.  20 

Sunnyvale,  California 

Hq.  AFCRL,  OAR  (CRRKP,  Ming  S.  Wong)  19 

L.  G.  Hanscom  Field,  Bedford,  Mass. 

G  68  Scientific  and  Technical  Information  Facility  1 

Attn:  NASA  Representative  (S-AK/DL) 

P.  O.  Box  5700 
Bethesda,  Maryland 


D-3 


I 

I 

I 

I 

I 

I 

I 


I 

1 

I 

1 

1 

I 


